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Abstract 

A G-strand is a map g{t, s) : M x M — )• G for a Lie group G that follows from Hamilton's 
principle for a certain class of G-invariant Lagrangians. The S'0(3)-strand is the G-strand 
version of the rigid body equation and it may be regarded physically as a continuous spin 
chain. Here, S'0(3)i^-strand dynamics for ellipsoidal rotations is derived as an Euler-Poincare 
system for a certain class of variations and recast as a Lie-Poisson system for coadjoint flow 
with the same Hamiltonian structure as for a perfect complex fluid. For a special Hamiltonian, 
the S'0(3)/^-strand is mapped into a completely integrable generalization of the classical 
chiral model for the ^©(Sj-strand. Analogous results are obtained for the S'p(2)-strand. The 
S'p(2)-strand is the G-strand version of the Sp(2) Bloch-Iserles ordinary differential equation, 
whose solutions exhibit dynamical sorting. Numerical solutions show nonlinear interactions 
of coherent wave-like solutions in both cases. Diff (M)-strand equations on the diffeomorphism 
group G = Diff (M) are also introduced and shown to admit solutions with singular support 
(e.g., peakons). 
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1 Euler-Poincare equations for a G-strand 

The Euler-Poincare (EP) theory of G-strands is an extension of the classical chiral models. The 
classical chiral models are nonlinear relativistically invariant Lagrangian field theories on group 
manifolds. As such, they are fundamental in theoretical physics. The vast literature of results for 
these models is fascinating. For example, it is well known that these models are integrable in 1 + 1 
dimensions and possess soliton solutions. See |5], [27] [351 iSSl [6], [261 1231 [2l] and references therein for 
discussions of the many aspects of integrability of the chiral models, including the famous dressing 
method for explicitly deriving the soliton solutions of these models, which is given in The 
solitons for the 0(3) chiral model are particularly familiar, because this model allows an integrable 
reduction to the one-component sine-Gordon equation (e.g. see [26]). Many generalizations of 
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these models have been introduced. For example, generalized chiral models with metrics that are 
not ad-invariant on the Lie algebra are considered in [29]. Other generalizations of chiral models 
for non-semisimple groups are studied in [16]. An integrable chiral model in 2+1 dimensions was 
proposed in [32] . Finally, the possibility for fermionic interpretatation of the current variables was 
explored in [33]. 

Left G-Invariant Lagrangian. We begin with the following ingredients of EP theory. For more 
details and in-depth discussion, see [IHIEI]- 

• Let G be a Lie group. A map g{t, s) : M x M — )• G has two types of tangent vectors, 
g := gt& Ug and g' := gs E Vg 

• Assume that the function L{g, Ug, Vg) : TG x TG — > M is left G-invariant. 

• Left G-invariance of L permits us to define / : g x g — )■ M by 



Conversely, this relation defines for any / : g x g — )• M a left G-invariant function L : 
TGxTG^ R. 

• For a map g{t, s) : M x M — G let 

X(t, s) := g'^gtit, s) = g~^g{t, s) and Y(t, s) := g'^gsit, s) = g~^g{t, s). 
Lemma 1.1. The left-trivialized tangent vectors X{t, s) and Y{t, s) at the identity of G satisfy 



Proof. The proof is standard and follows from equality of cross derivatives gts = gst, cf. [HI |2T] . 



KQ Ug,g Vg) = L{g,Ug,Vg). 



adxY. 



(1.1) 



As a consequence, equation (1.1) is often called a zero- curvature relation. 



□ 



Theorem 1.2 (Euler-Poincare theorem). 

With the preceding notation, the following statements are equivalent: 



i Hamilton's variational principle on TG x TG 




(1.2) 



holds, for variations Sg(t, s) of g(t, s 



) vanishing at the endpoints in t and s. 



ii The function g(t, s) satisfies Euler-Lagrange equations for L on G, given by 



dL d dL d dL 

dg dt dgt ds dgs 







(1.3) 
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iii The constrained variational principl^ 

S [ ' l{X{t, s),Y{t, s)) dsdt = (1.4) 
Jti 

holds on Q X g, using variations ofX and Y of the forms 

6X = Z + adxZ and 5Y = Z' + adyZ , (1.5) 
where Z{t,s) G Q vanishes at the endpoints. 
iv The Euler-Poincare equations hold on g* x g* 

d 61 _ 61 d 61 _ 61 ^ ,^ ^, 

dt 6X ^6X ds6Y ^ 6Y ^ ^ 

Proof. The equivalence of i and ii holds for any configuration manifold and so, in particular, it 
holds in this case. 

Next we show the equivalence of iii and iv. Indeed, using the definitions, integrating by parts, 
and taking into account that Z vanishes at the endpoints (ti, Si,t2, S2), allows one to compute the 
variation of the integral as 



6 



/(X(t, s), Y(t, s)) ds dt = (^^^sx'j + (^^,6Y^ ds dt 




/ 61 ■ \ / 61 \ 

— ,Z + adxZj + (— ,Z' + adyZ) dsdt (1.7) 

I d 61 ^,61 d 61 61 \ ^ 

( — — — ady -7^, + -TTz-, — ady ttt , Z ) dsdt . 
HiJsi \dt6X ^ 6X ds6Y ^ 6Y / 

To understand the notation in the last step, recall that the coadjoint action (ad* : x g* — )■ g*) 
is defined as the dual of the adjoint action (ad : x g — > g) with respect to the pairing 
(■,-):0*X0^Mby 



Thus, the calculation in (1.7) allows us to conclude that iii and iv are equivalent. 

Finally we show that i and iii are equivalent. First note that the left G-invariance of L : 
TG X TG — )■ M and the definitions of X{t,s) and Y{t,s) imply that the integrands in (1.2) and 
(1.4) are equal. Moreover, all variations 6g{t, s) G TG x TG of g{t, s) with fixed endpoints induce 
and are induced by variations {6X{t, s), 6Y{t, s)) G x g of X{t, s) and Y(t, s) of the form 

5X = Z + adxZ and 5Y = Z' + adYZ, (1.9) 

with Z(t,s) G g vanishing at the endpoints. The relation between 6g(t,s) and Z(t,s) is given by 
Z{t,s)=git,s)-^6g{t,s). 

Thus, if i holds, we define Z{t,s) = g~'^6g{t,s) for a variation 6g{t,s) with fixed endpoints. 
Then if we let X = g~^g{t,s) and Y = g~^g'{t,s), we have equation (1.9) by the same standard 

^As with the basic Euler-Poincare equations [25], this is not strictly a variational principle in the same sense as 
the standard Hamilton's principle. It is more like the Lagrange d'Alembert principle, because we impose the stated 
constraints on the variations allowed [H]. 
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calculation that produced (1.1). Conversely, if equation (1.9) holds with Z{t,s) vanishing at the 
endpoints, we define Sg(t, s) = g(t, s) o Z(t, s) and the above proposition guarantees then that 
this 6g{t, s) is the general variation of g{t, s) vanishing at the endpoints. Hence, i and iii are 
equivalent. □ 

Remark 1.3. For a right G-invariant Lagrangian, the results and proofs are the same as above, 
except for the sign change appearing in the variations, for which 6X = Z— adxZ and 6Y = Z adyZ. 



Evolutionary G-strand. We now define the fundamental quantity of interest in the remainder 
of the paper. 

Definition 1.4. A G-strand is an evolutionary map into a Lie group G, g{t,s) : M x M — )■ G, 
whose dynamics in t and s may he obtained from Hamilton's principle for a G-invariant reduced 
Lagrangian Z : g x g — )■ M, where q is the Lie algebra of the group G. The G-strand system of 
evolutionary partial differential equations for a left G-invariant reduced Lagrangian consists of the 



the zero-curvature equation (1.1) and the Euler-Poincare (EP) variational equations (1.6). 



Remark 1.5. Subclasses of the G-strand maps contain the harmonic maps into Lie groups studied 
mathematically in f31f and the principal chiral models of field theory in theoretical physics, reviewed, 
e.g., in [24^ . An interpretation of the G-strand equations as the dynamics of a continuous spin 
chain is given in fT^ . This is the source of the term, 'strand'. See also JE/- 



The evolutionary G-strand maps that we consider here arise from Hamilton's principle ^S* = 
with S = J £dt for a left G-invariant Lagrangian £ : g x g — )■ M given by 

pb POO 

S= i{X,Y)dsdt, (1.10) 

J aJ — oo 

with (X, Y) G X where q is the left Lie algebra of the group G. That is, X and Y are left-trivialized 
tangent vectors at the identity of the group, expressible as 

X{t,s) = g-^dtg{t,s) and Y{t, s) = g-^dsg{t, s) . (1.11) 

Remark 1.6. The distinction between between the maps (X, Y) : M x M — > g x g and their pointwise 
values (X(t, s), Y(t, s)) G g x g will always be clear in context, so that no confusion will arise. 
Likewise, for the variational derivatives 4^ and 4^ . 



Overview and organization of the paper. After this introduction, we shall begin developing 
G-strand dynamics in the example of the rotation group S0{3) in Section |2} In a certain case, 
this example recovers the well-known completely integrable principal chiral model [351 ES]- We 
extend these considerations in Section [3] to the isotropy group 5*0(3)^^ of a quadratic form defined 
by a symmetric matrix, K = K-^. In Section |4| the S'0(3)x-strand equations for a certain choice 
of Lagrangian are identified with the integrable dynamics associated with another variant of the 
principal chiral model, the completely integrable P-chiral model of [3l |3l] . Section |5] considers G- 
strand dynamics for Sp{2), the two-dimensional symplectic group, and extends the Bloch-Iserles 
ordinary differential equations to the case of partial differential equations for the S'p(2)-strand. 
By using the 'gauge' transformation that diagonalizes K and a series of canonical isomorphisms, 
the S'p(2)-strand equations for a particular choice of Lagrangian are also brought into the form of 
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the integrable S0{3)k P-chiral model. In each of these cases, we apply the Euler-Poincare (EP) 
framework for group invariant Lagrangians, then pass to the corresponding Lie-Poisson Hamil- 
tonian framework, where we are able to identify Hamiltonians that produce integrable dynamics. 
Numerical solutions for the 0(3) and Sp{2) versions of the P-Chiral models are provided in Section 
[6j In Section [7] we extend the G-strand framework to the case of diffeomorphisms on the real line 
and show that the Diff(]R)-strand dynamics admits singular solutions associated with a pair of 
momentum maps. Finally, Section [8] summarizes our results and provides some outlook for future 
research. 

In the standard Euler-Poincare (EP) framework, one applies the following steps [T8| [25]. 

(a) Write the auxiliary equation for the evolution of Y : M x M — )• g, obtained by differentiating 

its definition with respect to time and invoking equality of cross derivatives. 

(b) Use the Euler-Poincare theorem for left-invariant Lagrangians to obtain the equation of motion 

for the quantity di/dX : M x M — >■ g*, where q* is the dual Lie algebra. 

The resulting Euler-Poincare equation will be an evolutionary partial differential equation 
(PDE). We assume homogeneous boundary conditions on X{t, s), Y{t, s) and vanishing end- 
point conditions on the variation Z = g^^6g(t, s) G Q. 

(c) Legendre transform the Lagrangian to obtain the corresponding Hamiltonian. Differentiate 

the Hamiltonian and determine its partial derivatives. Write the Euler-Poincare equation 
in Lie-Poisson Hamiltonian form, in terms of the new "angular momentum" variable P = 
6i/SX G Q*. 

(d) Write the Lie-Poisson Hamiltonian formulation for G = S0{3) in terms of vectors. 

(e) Apply the EP procedure to other interesting choices of the Lie group G, e.g., 5^(2), and other 

choices of the Lagrangian i{X, Y) in Hamilton's principle. 

The EP framework with steps (a)-(e) provides the organization for each of the sections that follow. 

2 Euler-Poincare (EP) procedure for S0{3) chiral model 

This Euler-Poincare procedure in steps (a)-(e) above produces a series of results that are to be 
described in this section for the S0{3) chiral model. These are the following. 

(a) The auxiliary equation 

According to their definitions in X{t,s) = g^^dtgit, s) and Y(t, s) = g~^dsg{t, s) are Lie- 

algebra-valued functions over M x M. The evolution of Y is obtained from these definitions by 
taking the difference of the two equations for the partial derivatives (9(Y()f:,s) and dsX{t,s) while 
invoking equality of cross derivatives. Hence, Y evolves by the adjoint operation. 



dtY{t, s) - dsX(t, s) = YX - X Y = [Y, X] =: -adxY 



(2.1) 
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This is the auxihary equation we seek for Y{t,s). In differential geometry, this relation is called 
a zero- curvature equation, because it implies that the curvature vanishes for the Lie-algebra- 
valued connection one-form 

A = g'^dg{t, s) = Xdt + Yds . 



Equation (2.1 ) is also expressible in terms of the connection 1-form A as the Mauer-Cartan relation, 

dA + AAA = 0. (2.2) 
See, e.g., [7j for further discussion of such matters in differential geometry. 



When augmented by the divergence relation d * A = 0, namely, 

d*A = d{Xds + Ydt) = {dtX - 9,Y) dt A ds = , 



(2.3) 



equations (2.2) and (2.3) comprise the classical chiral model. Equation (2.3) may also be 
understood as a harmonic map from (t, s) G M x M into the Lie group G 



(b) The Euler-Poincare theorem 

According to the Euler-Poincare theorem, Hamilton's principle 5S = for S = J^^J^^ ^(X, Y) ds dt 
leads to the result in equation (1.7), 



d 



adv 



6e 



d 



dt6X ""^5X ds6Y^'''^^6Y- 
This is the Euler-Poincare equation for 6i/6X G Q* and it involves 6i/6Y G Q*. 



(2.4) 



Proposition 2.1. The Euler-Poincare equation (2.4) rnay he written as a conservation law for 

angular momentum P = 6£/6X; viz, 



d ( 51 



61 



(2.5) 



Proof. This formula follows from a standard result [T8j. Namely, let g{s) be a path with parameter 
s in a Lie group G and let Q(s) be a path in the dual q* of its Lie algebra q. Then the coadjoint 
operation Ad* : G x g* — )■ 0* satisfies 



Ad*,,._.Q(s) = Ad 



9{sr 



dQ 

ds 



adx(.)Q(s) 



(2.6) 



where X(s) = g{s) g{s) G Q. Hence, one finds the conservation law (2.5) by using (2.6) on either 



side of equation (2.4). 



□ 



Summary. In summary, the G-strand equations form a system of partial differential equations 



(PDE) comprising the Euler-Poincare equation (2.4) and its auxiliary zero-curvature equation 



(2.1 ), as follows. 



d 



+ 



d Si 



dt6X ds6Y 



adx^ + ad; 



6i 
dY 



(2.7) 



- adxY . 



(2.8) 
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Remark 2.2. The generalizations of the G-strand equations (2.1) and (2.8) to higher spatial 
dimensions are straightforward. 



Proposition 2.3. When the Lagrangian £(X, Y) : g — t- M given by 

i{X,Y) = l\X\'-l\Y\\ 



(2.9) 



then Hamilton's principle recovers the chiral model equation (2.5) 



Proof. The proof follows by comparing equations (2.4) and (2.3) for this choice of Lagrangian. □ 



(c) Lie-Poisson Hamiltonian formulation 

Legendre transforming the Lagrangian £(X, Y) : g x g — )■ M yields the following Hamiltonian 
h{P,Y) : Q* xg^R 

MP,Y) = (P,X\-£(X,Y). (2.10) 



Its partial derivatives are determined from 

6h 

' 5P,X) + (P 



Vanishing of the middle term in the second line defines the momentum P G g*. These derivatives 
allow one to rewrite the Euler-Poincare equation solely in terms of momentum P as 



^tP 

aY 



adx 



5h/SP 



p + a 



6h 



'6Y 

<9s^ — a-dsh/sp Y . 



— adt 



6h 
6Y 



(2.11) 



The corresponding evolution equation for any functional of /(P, Y) then follows as 



d_ 
di 



/(P,Y) 



ad 



>p + a 



5h 



— ad 



5f 



+ ( a^-ad 



Y 

Sh/5P^ , ^ 



6f 6h 
5P ' 6P 
5h 5f 
W 6P 



d. 



5f 5h 
6Y' 6P 



+ Y , ad^j/gp 
{/,M(P,Y) 



6h 



A* \ 

Sf/SP ^ - ^O^Sh/SP ^ / 
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in which the bracket { ■ , ■ } satisfies the properties defining a Poisson bracket because it is dual to 
the Lie bracket. Assembhng these equations into Hamihonian form gives, 



di 



'p' 




Y 





adgP ds — ady 
(9., + ady 



'6h/6P' 
5h/5Y 



(2.12) 



in which the box □ indicates how the ad- and ad*- operations occur in the matrix muhiphcation. 
In particular, 

ad^P(5V5P) = ad^V5pP, 
so each matrix entry acts on its corresponding vector component]^ 

Higher dimensions (G-branes) 



The notation in the Hamiltonian matrix (2.12) indicates how the jump to higher dimensions in 
the variable s G may be made naturally. This is done by using the gradient dj = d/ds\ 
J = 1, 2, . . . , to define the left invariant auxiliary variable Y,- = g^^djg in higher dimensions. The 



lower left entry of the matrix (2.12) defines a covariant gradient, and its upper right entry defines 
the adjoint operator, the covariant divergence. More explicitly, in terms of indices and partial 
differential operators, this Hamiltonian matrix becomes. 



d_ 
di 



Po 

Y? 



5. 



Sh/SP^' 
5h/5Yl 



(2.13) 



where the Hamiltonian structure matrix Baf3 is given explicitly as [T7] . 



B. 



a/3 



-P 



''/3k ' i 



A/35. I +/3 Y'^' 




(2.14) 



Here, the summation convention is enforced on repeated indices. Upper Greek indices refer to the 
Lie algebraic basis set, lower Greek indices refer to the dual basis and Latin indices refer to the 
spatial reference frame. The partial derivatives di = d/ds^, i = 1,2, ... ,n, act to the right on all 
terms in a product by the chain rule. 



(d) Write the Lie-Poisson Hamiltonian form for G = SO {3) 



For the case that t'^^ are structure constants for the Lie algebra so (3), then t^^^ = ea/sn with 
ei23 = +1 and the Lie bracket on so(3) may be identified with the vector product on M^. The 



Lie-Poisson Hamiltonian matrix in (2.14) may then be rewritten in vector form as 



d 


P 




Ft 







Px dj + YjX 

di + \iX 



'6h/6P' 
Sh/SYj 



(2.15) 



Returning to one dimension, one sees that stationary solutions for either dt — )■ 0, or (9^ — )■ satisfy 
equations of the same form as the heavy top. That the equations for either type of solution would 



^The matrix in equation (2.12) is also the lower right corner of the Hamiltonian matrix for a perfect complex 
fluid [H] . It also appears in the Lie-Poisson brackets for Yang-Mills fluids [12l [13] and for spin glasses [19] . See 
those references for full discussions of the Lie algebraic meaning of such Lie-Poisson brackets for the perfect complex 
fluid. 
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have the same form might be expected, because of the exchange symmetry under t ^ s and 
X f-)- Y. Perhaps less expected is that these equations would appear in the form of a heavy top. 
The EP formulation of the heavy top equations appears, for example, in [T8] . 



For G = S0{3) and Lagrangian £(X, Y) 



X 



—7- M, in one spatial dimension the Euler- 



Poincare equation and its Hamiltonian form are given in terms of vector operations in 
follows. First, the Euler-Poincare equation (2.4) becomes 



as 



di6X 



X X 



6£_ 
6X 



d_5i_ 



Y X 



6i_ 
5Y 



and its auxiliary equation (2.1) becomes 

dt 



d_ 
ds 



X + Y X X 



(2.16) 



(2.17) 



The Hamiltonian structures of these equations on so(3)* are then obtained from the Legendre 
transform relations 



5£ 
6X 



Sh 
6P 



X and 



6h 
6Y 



Hence, the Euler-Poincare equation (2.4) becomes 

d . 



dt 



P X 



6h 
6P 



d 6h 
d~s6Y 



and the auxiliary equation (2.1) becomes 

^Y. 

dt 



Y X 



Sh 



6i 
6Y 



6h 
6Y 



d 6h ^ 

I- Y X 

ds5P 5P 



(2.18) 



(2.19) 



which recovers the Lie-Poisson Hamiltonian form in equation (2.15). 

Finally, the equations for reconstructing 0(t,s) G S0{3) from the evolution of X(t, s) and 



Y(t, s) may be expressed by using the hat map 



so(3), e.g., Xij = -eijkX^ 



as 



dtO{t,s) = 0{t,s)X{t,s) and 
dsOit.s) = 0(t,s)Y(t,s). 



(2.20) 



Remark 2.4. The Euler-Poincare equations for the G-strand discussed here and their Lie-Poisson 
Hamiltonian formulation provide a framework for systematically investigating three-dimensional 
orientation dynamics along a fixed one-dimensional string. These partial differential equations are 
interesting in their own right and they have many possible applications. For an idea of where the 
applications of these equations could lead if the string were also allowed to bend self- consistently, 
consult fM W- 



The remainder of the paper will address the final item (e) in the Euler-Poincare G-strand 
procedure, by making two other choices for the group and the Lagrangian. 
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3 The S0{3)k strand 



Another interesting choice for the Lagrangian extends the orthogonal group 5*0 (3) to the isotropy 
group of a symmetric quadratic form in n dimensions. 

Proposition 3.1 (Isotropy group). Let K G GL{?),M.) he a symmetric matrix, K = . The set 

50(3)a' = {U e GL{3, R)\ U^KU = K} , (3.1) 
defines the subgroup S0{3)k o/GL(3,]R). Moreover, this isotropy group of K is a submanifold of 

Proposition 3.2. Let X = U^^U and Y = U^^U' for U ^ S, where U' and U denote derivatives 
ofU in s andt, respectively. Then the matrix quantities X := KX andY := KY are antisymmetric 

X^ = -X and Y^ = -Y. (3.2) 



Proof. This is a direct verification. 



□ 



Let X' and Y' denote separate derivatives of X and Y. Then, by an easy calculation one finds, 

X' - ?• = XY - YX =: [X,Y]. 
Hence, upon substituting the definitions of X and Y, one finds 

X' = A"X' = KY-+K[X,Y] 

= Y + Xi^-^Y - Yi^-^X 
=: Y+[X,Y]^ 
=: Y + adfY. 



(3.3) 



The corresponding auxiliary equation is to be compared with (2.1), 

dtY{t,s)-dsX{t,s) = -a.d^Y, 
in which Y evolves by the i^-adjoint operation. 

The same calculation of Hamilton's principle SS = a.s in the first section now yields the 



following Euler-Poincare equation, which is to be compared with (2.4), 

d Si 



ad 



K: 



Si d 



dtSX SX dsSY^ SY 



in which Si/SX evolves by the i^-coadjoint operation, defined as 

Si ^\ I 



adx , L 



\sx 



ad^Z 



(3.4) 



(3.5) 



Upon specializing to n = 3 and identifying the Lie bracket on so (3) with the vector product on 
as we did before, we find the corresponding /T- vector equations for the motion on the isotropy 
group S, by mapping 

adx Y -> ^(X X Y) and ad^*P ^ - X x i^P . (3.6) 
This means the Lie-Poisson Hamiltonian matrix in ( 2.14[ ) may be rewritten in iiT- vector form as 



d_ 
di 



KPxn ds + \x{Kn] 
a + K(\ X □) 



Sh/SP' 
Sh/SY 



(3.7) 



In the next section, the i^-vector equations obtained this way will be compared with the 
completely integrable P-chiral model [3, El] . 
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4 Relation to an integrable system: the P-chiral model 



For the Lagrangian £ = ^(|Xp — |Yp) in (2.9) the equations (3.3) - (3.6) in vector form are 



dtY{t, s) - dsX{t, s) + K{XxY) = 0, 



(4.1) 



dsYit, s) - dtX{t, s)+V X KV -Xx KX = 0. 

These may be cast into Lie-Poisson Hamiltonian form, as follows. 



(4.2) 



4.1 Hamiltonian formulation of the P-chiral model in 



X 



The Legendre transformation (2.10) in this case yields P = X and 

1, 



h{P,Y) 



2' ' 2' ' 



(4.3) 



which is sign-definite. Equation (3.7) now verifies the equations of motion in Lie-Poisson Hamil- 
tonian form, 



dtP = KP xP + ds\ + \ X {K\) , 
dtY = dsP + K(y X P) . 

We introduce a new diagonal matrix P, according to 

1 



P= ^ti{K)l-K , 



K = tT{P)l-P. 

With this definition one may easily check the identity 

K{X X Y) = X X PY - Y X PX . 



By using this identity, equations (4.1) - (4.2) can be written as 



dtY{t, s) - d,X{t, s) + X X PY - Y X PX = , 
d^it, s) - dtX{t, s) - Y X PY + X X PX = . 



(4.4) 
(4.5) 



(4.6) 



(4.7) 



(4.8) 
(4.9) 



4.2 Integrability of the P-Chiral model 

Under the linear change of variables 

u = X — Y and v = 



X- Y 



(4.10) 



equations (4.8) and (4.9) acquire the form of the following 0(3) P-Chiral model, [3], [31], [TT], 



dtu{t, s) + dsu{t, s) + u X Pv = - 



(4.11) 
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dtv{t, s) — dsy{t, s) + V X Pu = - 



(4.12) 



where the diagonal matrix P is defined in terms of the 3x3 symmetric matrix K in equation (4.6). 



The system (4.11) - (4.12) represents two cross-coupled equations for u and v. These equations 
preserve the magnitudes |up and |vp, so they allow the further assumption that the vector fields 



(u,v) take values on the product of unit spheres S x c 



The P-Chiral model is an 



integrable system and its Lax pair [5^ in terms of (u,v) is given in [3J. Expressing its Lax pair in 
terms of (X, Y) utilizes the following isomorphism between so(3) ©5o(3) and 5o(4): 



A(X,Y) 



/ 





Xs 


-X2 


Y^\ 









Xi 


Y2 











Y3 


V 




-Y2 




0/ 



(4.13) 



The system (4.1) - (4.2) can be recovered as a compatibility condition of the operators 



L = ds- A{Y,X){Xld + J), 
M = dt- A{X,Y){Xld + J), 



where the diagonal matrix J is defined by 
1 



J 



diag(-iri + K2 + Ks, K1-K2 + Ks, K1 + K2- K3, K1 + K2 + Ks 



(4.14) 
(4.15) 



(4.16) 



These steps get us to the 0(3) P-chiral model and explain its derivation as an Euler-Poincare 
equation and thus as a Lie-Poisson system. Previous derivations had identified a linear Poisson 
structure by other methods, mainly based on Lax pairs and Maurer-Cartan, or zero-curvature 
relations. Now that we have taken these steps, we can build a chiral-type model as an Euler- 
Poincare equation with quadratic kinetic and potential energy on any Lie group. 



Remark 4.1. // K = Id, equations {4-1) - {4-2) recover the 0(3) chiral model. 

5 The Sp{2) strand 

5.1 Bloch-Iserles (BI) equation 

The Bloch-Iserles equation is the n x n matrix differential system [2] 

X{t) = [N,X% 



(5.1) 



where the skew-symmetric matrix G so{n) is given and the solution X{t) G Sym(?T,) is a sym- 
metric matrix. The BI equation is expressible equivalently as 



X{t) = [B{X),X], with B{X) = NX + XN , 



(5.2) 



where B{X) : Sym(n) so{n). Consequently, the solution of (5.2) for X{t) is given by the 
similarity transformation 



X{t) = Q{t)X{(})Q-\t) , with QQ-\t) = B{X{t)) G so(n) . 



(5.3) 
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This form of the solution reveals that the BI equation is isospectral, i.e., the eigenvalues of the 
symmetric matrix X{t) are preserved. Isospectrality leads to constants of motion, obtained from 
the commutator form of the BI equation with B{X) = NX + XN, 

X{t) = [NX + XN + XN'^,X + XN], with XeR. (5.4) 

Consequently, the matrix invariants of X{t) are conserved, 

tr(X + AX)'' = const , for k = 1,2, . . . ,n - 1 . (5.5) 

For more information about the further analysis of the n-component BI equation in general form, 
see [2] and references therein. In what follows, we introduce the simplest case of the BI equation, 
for the symplectic group Sp{2), then we treat its corresponding G-strand extension. 



5.2 The simplest BI equation 

The simplest nontrivial form of the BI equation may be derived as a particular Euler-Poincare 
equation on the symplectic group Sp{2). Euler-Poincare equations on the symplectic group Sp(2n) 
have been discussed earlier in p3l [TO] and follow a familiar pattern. 

Let the set of 2 x 2 matrices Mj with i = 1,2,3 satisfy the defining relation for the symplectic 
Lie group Sp{2), 

MiJMf = J with ~ '^^ index z = 1, 2, 3. (5.6) 

The corresponding elements of its Lie algebra = MiM^^ E sp(2) satisfy {Jmi)^ = Jnii for each 
i = 1,2,3. Thus, Xj = Jrrii satisfying Xf = Xj comprises a set of three symmetric 2x2 matrices. 
For definiteness, we may choose a basis given by 

Xi = Jnii = ^ , X2 = Jma =(^92)' " "^"^^ " (l o) ' ^^'^^ 

This basis corresponds to the vector of momentum maps given by quadratic phase-space functions 
^ = (IqP) IpPjQ ■ p)^ that are often used in geometric optics and in the theory of particle beam 
design. That is, for z = (q, p)^ one identifies 

^z^XiZ = |q|2 = Xi, ^z^XaZ = \p\^ = X^, ^z%z = q • p = X3 . (5.8) 

Lemma 5.1. For X = Jm, Y = Jn E sym{2) with m,n E 5p(2), the J-bracket 

[X, Y]j := XJY - YJX =: 2sym(XJY) =: ad^Y (5.9) 
satisfies the Jacobi identity. 

Proof. By a straightforward calculation using = — Id2x2 with the definitions of X and Y, one 
finds 

[X, Y] J := XJY - YJX = - J(mn - nm) = -J[m, n] . (5.10) 

The lemma then follows from the Jacobi identity for the symplectic Lie algebra and linearity in 
the definitions of X, Y G sym{2) in terms of m,n G sp(2). □ 
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5.3 The BI Sp{2) strand 

Lemma 5.2. Let X = JMtM'^ with time derivative Mt = dtM{t,s) and let Y = JM'M'^ with 
derivative M' = dsM{t, s). Then 

X' = Y + [X,Y]j, (5.11) 



for the J-bracket defined in equation (5.9). 



Proof. Let m = MtM~^ and n = M'M~^. Then equality of cross derivatives implies the relation 

m' — n = nm — mn =: —[m,n] . 

Hence, upon substituting the definitions of the symmetric matrices X and Y, and using r = - J, 
one finds 

X' = Jm' = Jh — J[m, n] 

= Y + [X,Y]j = Y + XJY- YJX 
= Y + 2sym(XJY) =: Y + ad^Y. 

□ 



Equation (5.11) provides the dynamics of the quantity Y(s,t). Lemma 5.1 and Lemma 5.2 



now 



allow us to prove the following theorem for the derivation of the EP equation for Sp{2). 



or 



Theorem 5.3. Hamilton's principle 6S = for S = J^^ i{X,Y) ds dt implies the EP equation fi 
the symplectic Lie algebra 5p{2) , 

d 6i ( ,^,di\ d M ( , ^ 

a«< + 2^>'"(^><§x) + §;^Y + 2^>'-(^^aY) ("2' 

Proof. By a direct calculation, integrating by parts, substituting and rearranging yields 
5^ = , 5X\ + , 5y\ dsdt 

b 



\5X ' / \6Y ' 

( — , 9tZ + ad3^z) + ^ — , + ad^Jz) dsdt 



- dt— + s^di*— , Z) + {- ds— + sid'y*^, , Z) dsdt 



d 5i J, 5i d 5i J, 5i 



Here the formulas for the variations 6X and 6Y are obtained from the usual equality of cross 
derivatives. One defines the J-coadjoint action (ad"'*) as the dual of the J-adjoint action (ad'^) 
with respect to the trace pairing ( ■ , ■ ) of symmetric matrices by 

. Z) - (I . ad^z) ^ (1 . 2sy,„(XJZ)) ^ ( - 2sym(jx|) . z) ^ (5.13) 

□ 



DD Holm, RI Ivanov & JR Percival 



G-Strands 



16 



Remark 5.4. Consider the evolution equation (5.12) in the case that 



^(X,Y) = ^ / tr(X^)-tr(Y^) ds 



(5.14) 



where tr denotes trace of a matrix. (This is the motion equation for the Sp{2) chain with respect 
to the trace pairing of symmetric matrices.) 

For this Lagrangian, we have di/dX = X and di/dY = —Y, so the Euler-Poincare equation 
(5.12) becomes 



X + 2sym ( JX^) = Y' + 2sym ( JY^) , 
or, equivalently, in terms of the ordinary commutator 

X + [X, X J + JX] = Y' + [Y, Y.J + JY] . 



(5.15) 



(5.16) 



Remark 5.5. When Y is absent, equation (5.16) for the BI strand recovers the BI ordinary 
differential equation (5.2) for the case that N = J. 



Remark 5.6. Either of the equivalent forms (5.15) or (5.16) of the EP(symp) motion equation 
must be closed by using the evolution equation for Y , obtained in equation (5.11) from the compat- 
ibility of cross derivatives. Namely, 



or, equivalently, 



Y = X' - 2sym (X JY) 



Y = X' - (XJY - Y JX) . 



(5.17) 
(5.18) 



Equations (5.16) and (5.18) may be compared against corresponding equations (4-2) and (4.1), 
respectively. In this comparison, these systems seem to be dual to each other, as symmetric is dual 
to antisymmetric. 



5.3.1 Lie-Poisson Hamiltonian form of Sp{2) strand dynamics 



The Hamiltonian form of the Euler-Poincare equation (5.16) and compatibility equation (5.11) 
with its associated Lie-Poisson bracket may be found as before from the Legendre transform 

r. di , dh 
The Lie-Poisson bracket may also be obtained in the usual way as 



with /i(P,Y) = tr(PX) -£(X,Y) 



d 


P" 




Ft 


Y 





J* 

ds + ad 



ad^*P 



dg — ad- 




J* 



Sh/SP' 
5h/5Y 



(5.19) 



This is the J-bracket version of the Lie-Poisson Hamiltonian structure in (2.12); so its Jacobi 
identity will follow from that of the J-bracket discussed earlier. 

Here, as in (5.20), one defines the J-coadjoint action (ad"^*) as the dual of the J-adjoint action 
(ad"^) with respect to the trace pairing ( ■ , ■ ) of symmetric matrices. Explicitly, this is 



dh 

P,2sym(— JZ 



2sym(j^P 



(5.20) 
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Remark 5.7. Legendre transforming the Lagrangian (5.14) that is a difference of squares leads to 
a Hamiltonian that is a sum of squares 



/i(P, ^) = \j tr(P^) + tr(Y2) ds , 
where, as before, tr denotes the trace of a matrix. 



(5.21) 



5.3.2 The isomorphism between Sp{2) and SL{2) strand equations 



The Lie algebras sp(2) and st(2) are isomorphic. The G-strand equations in terms of the variables 
n e s[(2) and m G st(2) are 



m — n' — [m^, m] + [n^, n] = 0, (EP equation) 
11 — m' + [n, m] = (compatibility), 



(5.22) 
(5.23) 



in which Jm and Jn are symmetric. The first of these is the st(2) EP equation with Lagrangian 

1 



tr(n"^n + m"' m)ds . 



(5.24) 



Equations (5.22) and (5.23) are related to the previously considered i^T-chiral model (that is, before 
specifying how K relates to P) with 



K 




(5.25) 



where however this K is not diagonal. The relation is based on the following linear invertible map 
between vectors X e and st(2) matrices: 



A(X) 



XI 

2 V2 

^2 ^ 

y/2 2 



GS[(2). 



(5.26) 



Consequently, one finds 



and 



[A(X),A(Y)] 
A(K-X) 



[A(X),A(Y)] 
A(X X Y) 



A(K-X X Y) 

A^(X) 



A^(X X Y) 

([A(X),A(Y)])^ 



Thus, the constraint equation for the K-model (4.1) in the s[(2) representation becomes 

A(Y) - A'(X) = -A(K . X X Y) = -[A(X), A(Y)]. 



(5.27) 
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Likewise, the EP equation for the i^'-model (4.2) becomes 

A'fY 



A(X) 



-A(Y X K-Y) + A(X X K-X) 

-[A(Y),A(KY)]^+[A(X),A(KX)]^ 
-[A(Y),A^(Y)]^+[A(X),A^(X)f 
-[A(Y),A^(Y)] + [A(X),A^(X)]. 



(5.2^ 



Now upon identifying m = A(X) and n = A(Y) one finds that the s[(2) model equations 



(5.22)-(5.23) above are also integrable, because they are isomorphic to (4.1) - (4.2). 



Furthermore, one can make use of the following property for the action of K given in (5.25): 
the vector K ■ X G corresponds to the matrix — KXK G so (3), where as usual X is the so (3) 
representation of the vector X. The equations of the sp(2) model according to the above s[(2) 



identification can be written in vector form (4.1) - (4.2). In the so(3) representation, this becomes 



X+ [KXK,X] = Y' + [KYK,Y] and Y = X' - K[Y, X]K . 

The symmetric matrix K may be written as K = OKO^, with diagonal matrix K = diag( 
and orthogonal matrix given by 








note = 



T 



0- 



(5.29) 

-1,1,1) 

(5.30) 



In terms of the new variables X = 0^X0, Y = O^YO, equations (5.29) keep their form 
X+ [KXK,X] = Y'+ [KYK,Y] and Y = X' - K[Y, X]K . 



(5.31) 



However, now the matrix K is diagonal and this case is integrable, since it reduces to the P- 
chiral model. Thus, with the 'gauge' transformation that diagonalizes K and a series of canonical 



isomorphisms, the sp(2) EP equations with Lagrangian (5.31) have also been brought into the 
form of the integrable so (3) P-chiral model. 



6 Numerical solutions for 0(3) and Sp{2) P-Chiral models 



The equations for the 0(3) P-Chiral model (4.11-4.12) can be rewritten in the form 



atu(t,o + u(t,OxPv(t,o = o. 



(6.1) 



dMt,r])+v{t,r])xPu{t,r]) = 0, (6.2) 

for new wave-tracking variables 

^ = s — t , T] = s + t . 

This formalism shows explicity the nature of the cross-coupling between the natural physical 
variables, u and v. At fixed ^ (i.e. along a characteristic of fixed slope, c = 1), u may vary only 
when it interacts with nonvanishing v. Similarly, at fixed rj (a characteristic of fixed slope c = — 1), 



DD Holm, RI Ivanov & JR Percival 



G-Strands 



19 



V varies only with nonvanishing u. Note also, as previously stated, the quantities |u|(.^) and |v|(?7) 
are fixed along their respective characteristics. This in turn implies constraints on the evolution 
of the original geometric varibles X = (u — v)/2 and — (v + u)/2. 

In the static s frame, compactly supported solutions of the G-strand models consist of trains 
of generally fixed travelling profiles in u and v along the line at speed c = ±1 and undergoing 
rotation when they collide, or equivalently matched symmetric and antisymmetic pairs in X,Y 
undergoing the same behaviour. 



The model (6.1 6.2) is particularly tractable for numerical simulation since the absence of 



explicit spatial derivatives removes the numerical dispersion typical in time integration of advective 



operators such as appear in (4.11 -4.12). The forcings due to of the nonlinear interaction terms can 
be modelled through standard time discretization techniques. Here we choose to use the implicit 
midpoint rule using cubic spline interpolation over a pointwise discretization of ^ and t]. 

Figure [T] shows waterfall plot comparisons of the three components of the numerical solution 
solution u and of v for the standard Chiral equations, with matrix 

K = L 



and hence, from (4.6) 



P = I. 

The equations are integrated on a periodic domain in s with initial conditions 

u{0, s) = (0, A{s - 0, 25), 0), for s E [-0.25, 0.75] 

v(0,s) = 0.75), 0,0), for s G [0.25,1.25] (6.3) 

for a Gaussian profile ^(s) = exp(— 256s^) 

This choice of profile for the experiment is effectively arbitrary and in fact any suitably smooth 
profile could be substituted. The observed solution behaviour is as expected, with rightward, c = 1, 
travelling waves in u and leftward, c = —1, waves in v. At points of collision the nonlinear terms 
cause the waves counterrotate around each other, primarily into the third component due to the 
value of u X V at t = 0. 

Following the methodology and isomorphisms of the previous section we may also apply the 



same numerical model and analysis for the BI strand equations (5.16) by specifying matrices 





r-1 














1 







^ 





1 ) 






and setting the vectors u, v to map to and from symmetric matrices X, Y through 

U = (mi,M2,U3) = - \=X-Y 
2 \ Us U1-U2 J 

V = t"! , f2, ^3 = - = —X — Y 



Figure 2 shows the results of an integration in the BI mode of the vector initial conditions (6.3). 
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Discussion. The numerical results exhibit fundamentally the same behaviour, except that the 
isotropy breaking in P introduces preferred axes to the rotation so that transfers between the 
second and third components of u and v are faster than motions involving the first component. 
Hence, in this case, we observe a rapid growth in in comparison to over the course each 
collision. It is in part this directionality in the rotations which allows the Bloch-Iserles ordinary 
differential equations (ODE) to act a sorting process. 

Given the sorting nature of the Bloch-Iserles ODE (when y = and X is independent of s), 
it is natural to question whether this behaviour is observed on the spatially dependent strand. 
To address this question, we now consider the evolution of the initial conditions for two forms of 
twisted helical strand, namely 

_ 1 / sin(27rs) + cos(27rs) \ y_l/^00 

^ ~ 2 V sin(27rs) - cos(27rs) J ' ^ ~ 2 V 

, sin(27rs) cos(27rs) \ 1 / 



cos(27rs) -sin(27rs) y' ' 2V0 
i.e. under the mapping to vector space, 

Xi = (sin(27rs), cos(27rs), 0), = 0, 

X2 = (0, sin(27rs), cos(27rs)), Y2 = 0, 
or, in the P-Chiral physical variables, 

Ui = (sin(27rs), cos(27rs), 0), Vi = (— sin(27rs), — cos(27rs), 0). 

U2 = (0, sin(27rs), cos(27rs)), V2 = (0, — sin(27rs), — cos(27rs)). 

These are two of the simplest (in the spectral sense) states which can exhibits spatial dependence. 
In u-v space the solutions would consist, in the absence of the other component of counter-rotation 
helices. Figures |4] and [5] show the variation of locus of the eigenvectors of the various matrices 
with time. In both cases it will be noted that the oscillatory behaviour caused by the wave motion 
dominates, in this sense motions on the strand retard the sorting action of the BI operator as the 
algorithm can only operate when characteristics in u and v originating from the same point in 
coincide. 
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Figure 1: Array of waterfall plots showing the t and s dependence and componentwise comparison 
of P-Chiral variables u and v for the numerical solution of the initial value problem (6.3) in a periodic 
domain of unit length. Here the components of u are plotted in solid black in rows down the page and 
the components of v in dashed blue in columns across the page. Snapshots of profiles in s are shown 
at various times t with the vertical offset proportional to t. In the absence of the signal in the other 
variable each profile would translate forever at unit wavespeed, rightwards for u, leftwards for v. This 
motion would be independent of the actual profile chosen. At points of interaction the variables are 
rotated antisymmetrically, introducing a clear signal in the cross-component direction, u^, as well 
as a much less obvious response in ui and V2- The magnitudes \u\ and \v\ are conserved in frames 
travelling with the wave, even through the collision. 
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Figure 2: Array of waterfall plots showing the t and s dependence and componentwise comparison of 
Bloch-lserles variables u and v for the numerical solution of the initial value problem (6.3) in a periodic 
domain of unit length. Here the components of u are plotted in solid black in rows down the page 
and the components of v in dashed blue in columns across the page. Snapshots of profiles in s are 
shown at various times t with the vertical offset proportional to t. In the absence of the signal in the 
other variable each profile would translate forever at unit wave speed, rightwards for u, leftwards for v. 
Compared with the previous figure, the much larger signal in U3 post collision is obvious. This is due 
to the anisotropic nature of the Bloch-lserles strand operator. The magnitudes \u\ and \v\ are again 
conserved in each travelling frame. 
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Figure 3: Isometric plots of the helical initial conditions used in the Bloch-lserles strand problem shown 
in terms of the Ui,Vi variables. In the first panel one period of the spatial domain, s is mapped to 
the vertical axis, with the M'^ vectors Ui(s) and Vi(s) plotted at 25 locations (dotted blue and green 
lines respectively), linking the point s with points s + Ui(s) and s + Vi(s). The locus of these points 
is indicated in the solid lines. The projections of these loci onto the s-y, s-z, and y-z planes are also 
indicated in dashed black for u and dotted black for v. The initial condition for U2, V2 has fundamentally 
the same form, oriented albeit rotated to twist around the s-axis. Animations of the evolution of these 
visualizations are available as supplementary materials for numerical experiments presented in this paper. 
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1,0 -10 1.0 -10 

(a) (b) 

Figure 4: Isometric plots of surfaces generated by the eigenvalues of the symmetric matrices X{s,t) 
(Panel a) and Y{s,t) (Panel b) under the strand Bloch-lserles equations for helical initial conditions, 

u = (sin(27rs), cos(27rs), 0), v = (sin(27rs), cos(27rs), 0), s G [0, 1], 

on a periodic domain. Time dependence is plotted along the vertical axis, with eigenvectors 
Ai,A2 on the horizontal axes. Linear oscillatory motions dominate, with the locus of eigenvectors lying 
close to a circle, whose radius varies sinusoidally with time in both variables, with a clear half period 
phase difference. Animations of the time evolution of these solutions are available as supplementary 
materials for all numerical experiments presented in this paper. 




(a) (b) 

Figure 5: Isometric plots of surfaces generated by the eigenvalues of the symmetric matrices X{s,t) 
(Panel a) and Y(s,t) (Panel b) under the Bloch-lserles strand equation for helical initial conditions. 
Time is shown vertically In all cases the linear oscillatory motion of period one dominates. Animations 
of the evolution of these solutions are available as supplementary materials for all numerical experiments 
presented in this paper. 
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7 The G-strand PDE for G = Diff (M) 

A particularly interesting G-strand system arises when we choose G = Diff (M) and the Lagrangian 
involves the Sobolev norm. This case is reminiscent of fluid dynamics and may be written 
naturally in terms of right-invariant tangent vectors {t,s,x) and 'j{t,s,x) deflned by 



dtg = vog and dgg = ^ o g , 



(7.1) 



ad^ — 



ad! 



^7 



(7.2) 



adi.7 . 



where the symbol o denotes composition of functions. In this right-invariant case, the G-strand 
PDE system in equations (2.1) and (2.4) for G = Diff(]R) with reduced Lagrangian ^(z/, 7) takes 
the following form, which generalizes to G = Diff(]R'^) in any number of spatial dimensions, 

d 6i d 51 
dt ds 57 

d'j du 

dt ds 

These equations form a subset of the equations studied in [HI HZJ dSl EO] for complex fluids. They 
also form a subset of the equations for molecular strands studied in [8]. The latter comparison 
further justifles the name G-strands for the systems being studied here. These equations may 
be derived from Hamilton's principle for an afiine Lie group action, under which the auxiliary 
equation for 7 may be interpreted as an advection law. This interpretation is discussed further in 
Pi- 

Upon setting m = 51 /5v and n 
maps M X M G 
PDEs in (t,s,x). 



5£/57, the right-invariant G-strand equations in (7.2) for 



Diff(]R) in one spatial dimension may be expressed as a system of two 1-1-2 



mt + ris 
It - 



- ad*m - ad!,ra = -{i>m)^ - mz/^ - (772)^ - 717^ 



ad^z/ 



(7.3) 



The corresponding Hamiltonian structure for these Diff(]R)-strand equations is obtained by Leg- 
endre transforming to 

/i(m, 7) = (m, u) — £(z/, 7) . 



One may then write the m-7 equations (7.3) in Lie-Poisson Hamiltonian form as 



d 


m 




dt 


J. 





— adtm ds + ad* 



d. — ad^ 







6h/6m = V 
bhjb^ = —n 



(7.4) 



in which j£(I 
the real line 



This is the Lie-Poisson bracket dual to the action of the semidirect-product Lie algebra 

= X(M)(s)A^(Dens)(M) ©G((9,) 

is the space of vector flelds and A^(Dens)(]R) is the space of 1-form densities on 
plus a generalized 2-cocycle G{d^. In the 7 — z/ entry of the Hamiltonian matrix 
in (7.4), one recognizes the vector-field covariant derivative in s, and finds its negative adjoint 
operator in the z/ — 7 entry. For more details about how such Lie-Poisson Hamiltonian structures 
arise in complex fluids with flnite-dimensional order-parameter groups (broken symmetries) and 
full discussions of their Lie algebraic properties, see [HI [TTl [TH |30] and references therein. The 
differences in sign between the Hamiltonian matrix in ( |7.4 ) and those in (2.12) and (5.19) are due 
to the differences between right-invariance here and left-invariance in the other two cases. 
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7.1 Semi-stationary solutions 



Remark 7.1. For s -independent solutions, the Diff(]R) equations (7.3) reduce to EP equations 
for a Lagrangzan defined on T(Diff(M) x X(M))/Diffxo(M)), ^n which DiSxoW C Diff(M) zs the 
isotropy subgroup of the vector field parameter Xo{M.). Namely, for solutions that are functions of 
t and X, 



mt = - adlm + ad* (-n) = -{ym)^ - mv^ - (772)^. - 717^. 
7t = - ad^z/ = -z/7^ + 7z/^ . 



(7.5) 



The EP equations for the reduced Lagrangian £(z/, 7) represent ideal fluid dynamics with an advected 



vector field, 7. It's Lie-Poisson form is given by the Hamiltonian matrix in (7.4), without the 
generalized 2-cocycle C{ds). 



The semi- stationary equations (7.5) may be interpreted as a one- dimensional analogue of ideal 



incompressible magnetohydrodynamics (HMHD) in three dimensions /7| l21^. To compare the semi- 



stationary versions of (7.3) with the HMHD EP equations, we recall the latter in three dimensions 
for fluid velocity u G A^(]R^) and magnetic field B G j£(]R'^), which may be rewritten in the present 
geometric notation as 



Ut = — ad*u + ad^i? , 
Bt = —adsu , 



(7.6) 



with divu = and divi? = 0. One sees the similarity with (7.5) immediately. These HMHD 
equations arise as EP equations from Hamilton's principle with Lagrangian J2T] 



B) = IJ {u'-B^)d: 



X 



\B\ 



L2 



(7.7) 



Equations ( 7.6) may be written in Lie-Poisson Hamiltonian form ( 7.4) as 



d 


u 




— ad^u ad^ 




5h/5u = u 


dt 


B 




— ade 




5h/5B = B 



(7. 



with MHD Hamiltonian h{u, B) given by the sum of the kinetic and magnetic energies, 

h{u, B) = ^j{u^ + B^) d^x . 



(7.9) 



For further discussion of the EP theory underlying the equations for ideal fluid dynamics with 
advected quantities, see f21^ . 

7.2 Relation to the Camassa-Holm equation 



An interesting subcase of the system of semi-stationary Diff (]R)-strand equations (7.5) arises when 

■ in t( 



one selects the Lagrangian expressed only in terms of u as its H^ norm on the real line 

1„ „o 



m 5 



(7.10) 
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with vanishing boundary conditions, as |a;| — ?■ oo. In this case, m = u — Uxx, so the semi-stationary 
G-strand equations in (7.5) provide an extension of the completely integrable Camassa-Holm (CH) 
equation jl]. 

Si 

rrit = — ad*m = — (z/m)^ — with m = — = u — Vxx ■ ('''■11) 

Of 

Specifically, these modified G-strand equations reduce in the absence of 7-dependence to the CH 
equation, which admits singular solutions known as peakons |1] 

m{t, x) = J2 Maimx - Q\t)) , (7.12) 

a 

where we sum in a G Z over the integers, or over any subset of the integers. The peakon solution 



(7.12) of the CH equation may be understood as a singular momentum map, as discussed in 



7.3 Peakon solutions of the Diff (M)-strand equations 

With the choice of Lagrangian using the norm. 



{^.l) = \Mm-\\\l\\m. (7.13) 



the Diff(M)-strand equations (7.3) admit peakon solutions in both momenta m and n, with con- 
tinuous velocities v and 7. We state this result in the following theorem. 

Theorem 7.2. The Diff(]R) -s^ranc? equations {7.3^ admit singular solutions 

m(s, t,x) = Y^ Ma{s, t)6{x - Q"(s, t)) , n{s, t,x) = J2 ^a(s, t)6{x - Q\s, t)) , (7.14) 

a a 

z/(s, t,x) = K*m = Y^ Ma{s, t)K{x, Q") , 7(5, t,x) = K*n = J2 ^a(s, t)K{x, Q") , 

a a 

that are peakons for K(x, y) = le~'^~^' . These singular solutions follow from Hamilton's principle 
SS = for the constrained action S = J L(z/, 7, Q) dt given by 

S = j l(y, 7) + 5^ M,(s, t) {dtQ\s, t) - v{Q\ s, t)) + N,{s, t) {dsQ\s, t) - 7(Q^ s, t))ds dt . 



and they are peakons for the Lagrangian i'(z/, 7) given in equation (7.13). 
Proof. After replacing 

.iQ', s,t)^l - CTit, s)) aud ,(Q". . t) ^ . - Q'i, .)) . 

one computes 6S in Hamilton's principle 6S = for the Lagrangian £(z/, 7) as, 

+ [Y. ^) i^tQ^i^^ - ^^Q"^ ^' ^)) + H ^^-{s, t) {dsQ\s, t) - 7(g^ s, t)) ds dt 

a a 
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Hence, stationarity of the constrained action 5* in the statement of the theorem imphes 
Sp: ^ = Y^Ma{t,s)6{x-Q^{t,s)), 



6u 



^7 



SMa-. dtQ\s,t) = u{Q\s,t), 



(7.15) 



The first two equations recover the forms of the singular solutions from equation (7.14) in the 
statement of the theorem. These expressions are singular momentum maps, as discussed in detail 
in ||20j. The corresponding velocities are expressed as 



u{x, s,t) = K *m = J2 Mb{s, t)K{x, Q^) , 7(x, s,t) = K*n = Y^ Nb{s, t)K{x, Q^) , (7.16) 

b b 

for a symmetric, positive-definite kernel K{x,y). 



Inserting the forms of the solutions for the tangent vectors, or velocities u and 7 in (|7.16|) into 

dK{x, 



the compatibility condition in (|7.3|) yields 



dtNtis, t) + dsMb + ^(iV^M, - MfeiV, 



dx 



+ E K'' ^^^ll^'\ NbM^ - M,N^) = . 



(7.17) 



b,c 



dx 



Evaluating this formula at x = Q°'{s,t) then yields a relation among the solution parameters as 



(7.18) 



where we have introduced the matrix notation K""^ := K{Q"',Q^). The symmetric matrix K"-'' = 
K^"' follows from the choice of kernel and is assumed to be invertible. Invertibility of K"-^ in the 
compatibility equation (7.18) implies a relation 



dsM, - + Y,iNbM, - M,N,) — {K''' - K''){K-^),, = 



(7.19) 



which may also be regarded as an evolution equation for the Lagrange multiplier Ne{s,t). 

Equation (7.18), or equivalently equation (7.19) also arises from applying equahty of cross 
derivatives d'^^Q = d^^Q to the two constraint equations in (7.15) imposed by the Lagrange multi- 
pliers Ma and A^^^. This may be seen by substituting the two constraint relations into the equations 
for the cross derivatives, 

dlQ^ = s,t) + Y^ ^^7(Q^) = daiQ", s,t) + Y, ^^KQ^) = dlQ , (7.20) 
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whereupon rearranging then yields equation (7.18) and invertibihty of K""^ imphes (7.19). 



Finally, the SQ" equation in the set (7.15) is given by 



5, + J] M, — + 9. iV. + iV. J] iV, — = . 



(7.21) 



which is an evolution equation for the Lagrange multiplier Ma{s,t). 



For the Lagrangian in (7.13) we have u = K * m and = K * n for the symmetric positive 
definite kernel K{x,y) = |e~'^~^', which means that in this case matrix elements of K^'' take the 
peakon shape K'^^ = ^e"''^""*^'''. The corresponding singular solutions (7.14) are peakons for the 
Lagrangian in (7.13). □ 



Summary. Upon collecting equations, we have the following three results. 

1. The singular solutions of the G-strand equations for G = Diff(]R) are represented by two 
singular momentum maps. 



m 



n 



s,t,x) = ^M„(s,t)(5(x-Q"(s,t)), T;Emb(§,M) ^ X*( 

a 

s,t,x) = J2Na{s,t)6{x - Q^{s,t)) , T;Emb(§,M) ^X:( 



(7.22) 



where § G M is the support set of the delta functions, and Emb(§, M) denotes the set of 
smooth embeddings Q : E> ^ M.. 

2. The two tangent vectors, or velocities, that correspond to these momentum maps are, 

z/(t, s,x) = K *m = Y^ Ma{s, t)K{x, Q^is, t)) , 

a 

7(t, s,x) = K*n = J2 Nais, t)K{x, Q"(s, t)) . 



(7.23) 



These arise from the convolution K* : j£*(]R) — ?■ X{M.) of m and n with the kernel K. 
3. The solution parameters {Q'^{s,t), Ma{s,t), Na{s,t)} with a G Z that specify the singular 



solutions (7.14) are determined by the following set of evolutionary PDEs in s and t, in 
which we denote K^^ := K{Q'', Q^): 



b 



dtMa{s, t) = - dsNa - 5^(M,M, + NaN,)^^ (no sum on a). 



(7.24) 



dtNeis, t) = dsMe + Y,{N,M, - MbN,) — {K^'' - K''){K-')ea (do sum on a). 



a,b,c 
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The last pair of equations in (7.24) may be solved as a system for the momenta, or Lagrange 
multipliers (M^, Ne), then used in the previous pair to update the support set of positions s). 
Given Q"(t, s) for a G Z, one constructs {m,n{t, s, x)) along the solution paths x = Q"'{t, s) from 
equations ( |7.22 ) and then obtains (z/, 7(t, s, x)) for x G M from equations (7.23). Alternatively, 
knowing the position Q°'{s,t), a = 1, . . . , n, for all s at a given time t, also determines Na upon 
inverting the matrix K""^ in the second equation in (7.24). 



7.4 Canonical Hamiltonian form of the Diff (M)-strand peakon dynamics 



The Diff(]R)-strand peakon equations (7.24) may be written in a canonical Hamiltonian form, after 
making a Legendre transformation to a Routhian, or constrained Hamiltonian, by forming 



H{M,Q) 



J2MadtQ^]ds-L{u,j,Q) 

^ Mau{x)6{x - Q'')dx - 7) - 5^ NaidsQ"" - 7(Q")) 



(7.25) 



ds 



where L(z/, 7,(5) denotes the constrained Lagrangian appearing in Theorem 7.2 
L{u,7,Q) ■■= [ ku,^) + J2{Ma{s,t),dtQ%s,t)-u{Q'^,s,t)) 



(7.26) 



ds . 



The variations of the resulting Routhian H are given by 

6i 



SH{M,Q) = j (Y.Ma5{x-Q% 



5u 



6i 



J2nJ{x-Q'^)~-,6^ 



(7.27) 



a, 6 



ds , 



where the angle brackets ( ■ , ■ ) denote pairing in x G M. Vanishing of the first two terms recovers 
the two singular moment maps in (7.22). Vanishing of the third term imposes the second constraint 
in (7.24). Finally, the remaining two terms produce the corresponding canonical equations for the 
Routhian H with canonically conjugate variables Q°'{s,t) and Ma{s,t), 



dtMa 



6Ma 

SH 



(7.28) 



Compatibility of the constraint dsQ"^ = 7(<5") and the canonical equation dtQ"" = viQ"") then 
recovers the last equation in (7.24) for the evolution of the Lagrange multiplier A^^- This completes 
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the canonical Hamiltonian formulation of the singular solutions of the Diff(M)-strand equations in 

(7.29) 



(7.24). The conserved energy of the peakons for the Lagrangian in (7.13) is 



E 



J2 MaK'^^Mh + NaK^'^Nb ds . 

a,b 



7.5 Alternative forms of the Diff (M)-strand equations 



Velocities, Eulerian. The Diff (]R)-strand equations (7.3) may be written equivalently in terms 
of the velocities v and 7, in Eulerian form as 



dtf + ffx + <9s7 + llx = - dxP{i^, 7) , 

where P(z/, := K* (z/^ + + j . 



(7.30) 



Velocities, Lagrangian. Equations (7.30) themselves may be rewritten equivalently by using 



Lagrangian time derivatives in s and t, in the compact form, 

= 0. 



du 


d'~) 


dt 


d'y 


dv 


'dt 


ds 



(7.31) 



Momenta, Lagrangian. Going back to the momentum 1-form densities m dx^ and n dx"^ allows 
us to rewrite equations (7.3) once more equally compactly using Lagrangian time derivatives as 

d 



dt 



(mdx ) + 



dx/dt=v 



d 

ds 



d 
dt 



{n dx^) 

dx/ds='y 
d 







7 - 

dx/dt=v 



ds 



(7.32) 



V = 

dx/ds=^ 



Outlook: Studies of the solution behavior and the issue of integrability of the G-strand equations 



(7.3) for G = Diff(]R) with the Lagrangian (7.13) may be expected in future investigations along 



the present lines. For example, the velocity form of the Diff (]R)-strand equations in (7.30) may 
be useful in determining whether they produce wave breaking from smooth initial conditions, by 
studying their slope dynamics at a moving inflection point of both velocities to prove a steepening 
lemma. 



Higher spatial dimensions. The G-strand equations may also be written in higher spatial 
dimensions. For example, the G-strand equations for volume-preserving diffeomorphisms in three 
dimensions with G = SDiff(]R^) may be written in terms of divergence-free vector fields as 

dtUJ + [u,uj] + dsfl + ['y,^!] = , with u) = cmlv and 17 = curl 7, (733) 
dtiy — dg^ — [j/, 7] = , where div = and div 7 = 0. 
It may be interesting to study the point-vortex dynamics of these equations in two dimensions. 
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8 Conclusions 



This paper has provided a general framework for studying the G-strand equations, the PDE system 
comprising the Euler-Poincare (EP) variational equations (2.4) for a G-invariant Lagrangian, cou- 
pled to the auxiliary zero- curvature equation (2.1). The latter has often been the departure point 
and main focus in other approaches, because it sets up the Lax-pair formulation of the system. 

Once the G-invariant Lagrangian has been specified, the EP and zero-curvature system of G- 
strand equations (2.4) and (2.1) follow automatically. In the present paper, the well-known 5*0(3) 
chiral model PDE system has been extended by placing it into the G-strand framework and then 
deriving its Lie-Poisson Hamiltonian formulation on the dual Lie algebra so(3)*. The corresponding 
Lie-Poisson bracket given in equation (2.12) turns out to be the same as for a perfect complex fluid 
[T7] . It also appears in the Lie-Poisson brackets for Yang-Mills fluids [12l |T3] and for spin glasses 
|19] . The SO (3) chiral model itself emerges in the G-strand framework for the very simple choice 



of quadratic Hamiltonian given in equation (4.3). 

In this framework, two new integrable 1+1 EP PDE systems have been identified in Section [3] 
and Section [5] that may be transformed into the 5*0(3) P-chiral model. These are the 5*0(3)^: and 
Sp{2) G-strand equations in which the corresponding Hamiltonians are sums of squares. Other 
types of integrable P-chiral models are also available. The various Lie algebras for these may be 
classified according to the signature sgn(K) of the real symmetric matrix K, defined as the triple 
(/, m, n) of the numbers of its positive, null, and negative eigenvalues, respectively. 

Finally, the dynamics of the soliton solutions for these two integrable systems has been studied 
numerically and their coherent nonlinear scattering behavior has been explored for some simple 
initial conditions in Section [6l 

Other systems may be treated using the same approach. For example, one may formulate the 
G-strand equations for the two-dimensional Toda lattice, in terms of an appropriate Lie-algebra 
of lower-triangular matrices, as reviewed in [15]. The multicomponent Bloch-Iserles system may 
also be formulated in this way. These other systems are likely to be fruitful subjects of future 
investigations of the G-strand PDE, including the G-strands on the diffeomorphisms, introduced 
in Section [71 
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